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ABSTRACT 

The most massive star clusters include several generations of stars with a different chemical composition (mainly revealed by an Na-0 anti- 
correlation) while low-mass star clusters appear to be chemically homogeneous. We are investigating the chemical composition of several 
clusters with masses of a few 10 4 M to establish the lower mass limit for the multiple stellar population phenomenon. Using FLAMES® VLT 
spectra we determine abundances of Fe, O, Na, and several other elements (a, Fe-peak, and neutron-capture elements) in the old open cluster 
Berkeley 39. This is a massive open cluster: M ~ 10 4 M Q , approximately at the border between small globular clusters and large open 
clusters. Our sample size of about 30 stars is one of the largest studied for abundances in any open cluster to date, and will be useful to 
determine improved cluster parameters, such as age, distance, and reddening when coupled with precise, well-calibrated photometry. We find 
that Berkeley 39 is slightly metal-poor, <[Fe/H]> = -0.20, in agreement with previous studies of this cluster. More importantly, we do not 
detect any star-to-star variation in the abundances of Fe, O, and Na within quite stringent upper limits. The r.m.s. scatter is 0.04, 0.10, and 
0.05 dex for Fe, O, and Na, respectively. This small spread can be entirely explained by the noise in the spectra and by uncertainties in the 
atmospheric parameters. We conclude that Berkeley 39 is a single-population cluster. 
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1. Introduction 

Most, perhaps all, Galactic globular clusters (GCs) host mul- 
tiple stellar populations that can be identified b y variations of 
Na and O abundances (e.g, Carretta et al. . 2010l and references 
therein), which are anti-co rrelated with each other (see e.g., 
iGratton et al. . 2004 . 2012 , for recent reviews). On the other 
hand, Na-rich and O-p oor stars that are common in GCs are 
rare among field stars ( Carretta et al.[) and are possib ly absent 
from open clusters (OCs, see e.g. Jde Silva et al. , 2009). A sim- 
ilar result i s obtained for other lig h t elements, such as C an d 
N, see e.g jMartell & Smithl d2009l) . lMartell & Grebell ( 1201 Oh . 
The onset of the mechanism that leads to the multiple stellar 
population phenomenon is tied to the dynamics of the star for- 
mation process and has in turn to depend on many variables 
such as age, metallicity, total mass, and probably also on the en- 



vironment (e.g., position in the Galaxy or in a smaller galaxy ) 



Send offprint requests to: A. Bragaglia, angela.bragaglia@oabo.inaf.it 
* Based on observations collected at ESO telescopes under pro- 
gramme 386.B-0009 

** Tables 2 and 3 are only available in electronic form at the CDS 
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where the cluster formed. According to iCarretta et al . (201 
however, there seems to be a minimum (present-day) threshold 
mass of about a few 10 4 M , as shown in Fig.[T](a n updated ver- 



sion o f the My-age plot originally presented in ICarretta et al 
1201 ()l) . 

Unfortunately, few clusters with masses close to this thresh- 
old value have been surveyed extensively enough to show 
whe ther they really have Na and O abunda nce variations (see 
e.g., iCarretta et all 1201 fj iBragaglial 1201 2b . Two GCs in this 
mass interval, Ter 7 and Pal 1 (Fig. [TJ seem to lack these 
variations; howeve r, only a handful of stars in each of them 
has been s tudied ( Tautvaisie ne et al. L 12004 Sbordone et al 



200 5L 120071: ICoheni 120041). Another case seems to be Rup 106 
dGeisler & Villanovai 120121) . at a higher present-day mass. A 
possible case for variations in the abundances of these ele- 
ments has r ecently been made f or the massive old open cluster 
NGC 6791 (iGeisler et al.Ll2012l but see the conflicting conclu- 
sion reached by our group, Bragaglia et al, in preparation). 

Observation of a larger sample of clusters is required, with 
the same kind of large samples we obtained for more than 
20 GCs using the multi-object, high-resolution spectrograph 



2 



A. Bragaglia et al.: The open cluster Be 39 



M v 
- 1 



-8 



-2 



_ 2x1 6 M e 


i i i i i i 

* * ~ 


I_2xJ0 5 M s 


* — j&m^, 


_ Ter7 W 




NGC6791 

-o o o 
* o 


a ; 

2xl0 3 M~ 


[ I I I I I I I I I 1 


0.6 


0.8 1 
Relative age 



Fig. 1. Relative age parameter vs. absolute magnitude My for 
globular and old open clusters. Red filled symbols are GCs 
where the Na-O anti-correlation has been observed; open star 
symbols and open squares mark MW and Sagittarius GCs, re- 
spectively, for which not enoug h data are availabl e; open circles 
are old open clusters (data from lLata et al . 2002). Finally, light 
blue crosses indicate several clusters that do not show evidence 
of Na-0 anti-correlation: two GC members of Sagittarius dSph, 
one of the main body (Ter 7) and one of the stream (Pal 12); 
the new entry Rup 106; and two massive, old OCs, Be 39 and 
NGC 6791 (the latter indicated also in red, since its situation 
is still not completely assessed). Two orange, large circles in- 
dicate the other GCs of our project . Superimposed are lin es of 
constant mass (light solid lines, see lBellazzini et al.l l2008). 



FLAMES @VLT (see e.g.. ICarretta et all l2009albh . We have 



started systematic observations of clusters in this transition re- 
gion, including both the most massive metal-rich disc open 
clusters and the faint end of the mass distribution of metal-poor 
globular clusters. This approach may help to understand if the 
only factor at play is mass (remembering of course that My 
measures the present-day mass, not directly the original one), 
or if other factors, such as the environment in which the clusters 
formed, have to be considered. 

In this paper we present results for the OC Berkeley 39, 
an old (about 6 Gyr, iKassis et alJ 1997 and Sect. 12. U mas- 
sive OC (My = -4.28. lLata et al.l " 20021 indicating a mass 
of ~ 2 x 10 4 M , see Fig. []])■ Be 39 is located in the outer 
part of the disc, a t about 1 1 kpc from the Galactic centre (see 



Friel et all [2oTob . Friel et al. determined the chemical com- 



position of this cluster from high-resolution spectra of four 
red giants. These authors showed that this cluster is slightly 
metal-poor ([Fe/H]=-0.21 +0.01), which agrees very well with 
the metallicities of other clusters at similar Galactocentric dis- 
tances (see their Fig. 3). Friel et al. obtained very consistent 



results from their four stars. However, their sample is too small 
to exclude a significant spread in abundances in the cluster, as 
would be expected from the multiple population phenomenon. 
Extant photometric data are not sufficiently accurate to shed 
light on this question (see Sect. 12. U . Chemical inhomogeni- 
eties in Be 39 cannot be ruled out at present. Therefore we 
included Be 39 in our survey, and obtained abundances from 
high-resolution spectra for several tens of member stars, a sub- 
stantial fraction of the evolved population of this cluster. 

While the spectral regions of our data were optimised 
to discuss the Na-O anti-correlation, our spectra have a 
wider interest. Open clusters are good tracers of the prop - 



erties of the disc of o u r Galaxy, s ee e.g ., iFriell (11995b : 



Bragaglia & Tosil (2006); lYong et all (1201 2l) and references 



therein. Notwithstanding many observational efforts, only a 
fraction of the k nown Galactic OC s (about 2000, according to 
the catalogue in iDias et all 120021 and its updates) have been 
observed with high-resolution spectroscopy. Searching the lit- 
erature, we found only about 80 with an age older than 100 
Myr (see Sect. |4| and in general only a handful of stars were 
observed in each of them. If one aims to investigate the history 
of the Galactic disc (formation and evolution), the old clusters 
are very important; our observations provide further, improved 
information on one of the few very old OCs. The large sample 
of stars is also very useful, because it permits us to put stringent 
constraints to the internal abundance homogeneity, something 
that is relevant for the cluster formation mechanisms. 

The paper is organized as follows: in §2 we present the ba- 
sic data for this cluster, the selection of stars to be observed, 
and the observations themselves. In §3 we describe the analy- 
sis method; a critical point is the derivation of accurate atmo- 
spheric parameters. In §4 we present and discuss the results of 
our analysis. §5 contains the summary and conclusions. 

2. Observations 

2. 1 . Photometry and cluster parameters 

Photometric observations for Be 39 have been used to study its 
variable stars and to determine the cluster para meters (see e.g., 
iKaluznv & Richtlerl Il989l; IKassis et all Il997h : a summary of 
past results can be found in iFriel etal.l(l2010l) Of t he pub licly 
available optical photometry, IKaluznv & Richtlerl (I1989I) col- 
lected BjJ^data^n a sm all field of view (4.5x7.2 arcmin 2 ) 
while Kassis et al. (1997) provided B, V, I photometry on a 
14.7x14.7 arcmin 2 area, bett er suited to the FLA MES spectro- 
graph field (25' in diameter. lPasquini et alll2002h . We obtained 
the two photometric catalogues through the database WEBDA0 
and used the second set to select our targets (see Fig.|2|. 

In our series of papers on the Na-O anti-correlation in 
GCs we used atmospheric parameters directly derived from 
the photometry, so we checked if it was possible to do the 
same for Be 39 as well. Unfortunately, the two photometric 
data sets are not identical; they show an average difference of 
-0.02 (rms=0.07) mag in V and +0.06 (rms=0.10) mag in B 
(in the sense Kassis et al. minus Kaluzny & Richtler) based 



The OC database is at http://www.univie.ac.at/webda/webda.html 
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Fig. 2. Left: CMD for Be 39 jKassis et allll997h with observed stars indicated by different symbols: UVES target (blue circles), 
GIRAFFE targets (orange squares for members, cyan triangles non-members), binaries (magenta stars). Right: CMD for the inner 
3.5 arcmin, with observed member stars indicated by larger symbols and three isochrones shown. The parameters used for the fit 
are indicated. 



on about 730 stars in common. This could seem a minor ef- 
fect, but a difference of about 0.07-0.08 mag in B - V means 
a large difference in stellar effective temperatures in Be 39. 
Furthermore, the largest discrepancies between the two sets are 
found also among the brightest stars, especially the red giant 
branch (RGB) and red clump (RC) ones. We were unable to 
unequivocally establish which, if either, of the tw o photomet- 



ric set s is the best calibrated, so we adopted the Kas sis et al 



dl997l) one because it is available over a larger field of view. 
However, we manually adjusted a few values for stars that are 
radial-velocity (RV) members of the cluster but have magni- 
tudes and colours that are discrepant from the evolutionary se- 
quences (see Sect. 13. U in one of the CMDs. 

We also re-deter mined the cluster p arameters, using the 



2008) retrieved from the 



Padova isochrones dMarigo et all 
interactive database http://stev.oapd.inaf.it/cgi-bin/cmd at the 
metallicity derived by iFriel et al.1 d2010b . i.e., [Fe/H]=-0.20. 



Fig.|2](right panel) shows the result; the age is between 5.5 and 
7 Gyr, with a best-guess age of about 6 Gyr. The distance mod- 
ulus is (m-M) () = 12.94, and the reddening is E(B-V) = 0.17, 
slightly higher than found in previous studies (but consistent, 
considering the relatively large uncertainties in the photometric 
data). These values were used to determine effective tempera- 
tures and surface gravities of the programme stars (Sect. [XT). 
Note however from the figure that RGB stars do not define a 
narrow sequence; this could be partly because of undetected 
binaries, but could also indicate differential reddening and/or 
uncertain photometric values. 

We could not locate Stromgren colours for Be 39, which 
could have been used to test the uni- or multi-population 
nature of the cluste r, as discussed on empiri cal grounds by 
Carretta et al.l Q201 lb and on theoretical ones bv lSbordone et al 



(1201 II) . The latter paper also considered the more usually avail- 
able Johnson-Cousins bands; however, as we have noted, ex- 
isting broadband photometry in those filters is not sufficiently 
precise for this goal. 

2.2. FLAMES data 

We obtained six exposures of Be 39 with the multi-object spec- 
trograph FLAMES @VLT, as in our previous studies on the 
Na-O anti-correlations in clusters. The observations were ob- 
tained in service mode and were spaced in time (see Tab. [TJ 
to search for RV variations that would indicate the presence 
of binary systems. We used the UVES 580nm setup (AA 
4800 - 6800 A) and the GIRAFFE high-resolution grating 
HR13, which contains the [O i] line at 6300 A and the Na i 
doublet at 6154-6160 A. 

We sele cted Be 39 candid ates using avai lable information 



on RV fromFriel et al. (2002) and lFriel et al.1 (2010) to help se 



lecting only probable cluster members. All chosen targets are 
free from neighbours: they have no other stars closer than 3" 
(or 2" if they are at least 2 mag fainter than the desired can- 
didates). As shown in Fig. [2] we selected stars on the RC, or 
close to it, for the UVES fibres; only one configuration was 
prepared, so we have seven stars observed at high-resolution 
(R =: 45000), with the eighth fibre used for sky subtraction. The 
GIRAFFE fibres (at R 22500) were allocated to other RC and 
RGB objects and to stars on the subgiant branch (SGB) and the 
turn-off (TO) of the main sequence; 16 fibres were used for 
sky correction. For the SGB and TO stars the spectral region 
is not ideal, since the [O i] line is too weak to be measured. 
The Na i features of these stars still could be used to search for 
a possible Na abundance spread. However, the main targets of 
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our observations are RGB and RC stars, for which Na and O 
can be measured; the most interesting information provided by 
fainter stars is their RV (i.e., membership). Information on all 
observed sta r s can be found in Tab. ??; we give the ID in the 
Kassis et al. I dl997h catalogue, the ID in WEBDA, equatorial 
coordinates, V, B — V and V - I colours (with a flag indicating 
stars corrected as described above), and the heliocentric RV. 

The spectra were reduced (bias and flat field corrected, 1- 
D extracted, and wavelength-calibrated) by the ESO personnel. 
We applied sky subtraction and division by an observed early- 
type star (UVES) or a synthetic spectrum (computed at the 
GIRAFFE spectral resolution) to correct for telluric features 
near the [O i] line, using the IRAI0 routine telluric. The latter 
correction was applied only to the UVES and bright GIRAFFE 
samples (RC, RGB stars), since this [O i] line is not usable 
in the warmer faint GIRAFFE stars (SGB and TO). We mea- 
sured the individual RVs using midlines in IRAF and shifted 
all spectra to zero velocity before combining them. The UVES 
final spectra have signal-to-noise ratio (S/N) in the range 100- 
170. The S/N ratios (computed on small continuum regions) 
of GIRAFFE spectra of RC stars are in the range 150-210; of 
RGB stars are generally between 100 and 200, and those of 
SGB and TO stars are 25-70. Values for individual stars are 
listed in Tab. ??. A portion of the UVES spectra showing the 
Na i lines region is shown in Fig. [3] 

The cluster average RV was computed separately for UVES 
and GIRAFFE spectra. We found < RV UVES >= 57.35 
(rms=0.85) km s" 1 and < RV C1RAFFE >= 58.21 (rms=1.62) 
km s _1 . We assume that all stars within ±3<x from the aver- 
age are cluster members, with some more dubious cases indi- 
cated in Tab. ??. Bo th these va l ues co mpare well to RV=58.65 
(rms=2. 34) km s' 1 dFriel et al.Ll2010i) and RV=55.0 (rms=2.9) 
km s" 1 ( Frinchabov et all 2006 ). Fig. |4] shows the RV distri- 
bution of our stars and compares it to the expected distribu- 
tion of field star ve locities, according to the Besan^on model 
dRobin et al. ■ l2003h . ie 39 stands out conspicuously and fewer 
than about ten field stars have RVs similar to the cluster one. 

After pruning the sample using RVs we found a total of 67 
possible single member stars. All seven UVES stars are RV 
members; of the GIRAFFE stars, 60 are possible members, 
eight are binary systems, and 43 are non-members. This is one 
of the largest sample of spectroscopically studied stars in OCs. 

While measuring the RVs we found a few stars with dis- 
crepant values at different observation times, which is indica- 
tive of binary systems. On closer inspection, some of them 
also displayed changes in spectrum appearance; an example is 
shown in Fig. [5] Th ere is also a suggestion, combining our RVs 
with those of Friel et ail d2010h (see Sect. [4/2}, that one of the 
UVES and one of the GIRAFFE targets may be binary. If all 
these systems are members, this indicates a binary fraction of 
10/75 or 13%. This is only a lower limit and it perfectly con- 
forms to the fractio n of binaries usua l ly fou nd in OCs from 
photometry (see e.g. jBragaglia & Tosil 120061) . 




6152 6154 6156 6158 6160 6162 
Wavelength (A) 

Fig. 3. Portion of the UVES spectra in the Na i lines region; the 
two Na lines are indicated by the dotted lines. The normalised 
spectra are offset for clarity and are plotted in order of decreas- 
ing r e ff, with the RGB star at the bottom. The Na and all other 
lines are practically identical in all stars, with only a small dif- 
ference between the RGB and the RC ones. 



Table 1. Log of the observations 



2 IRAF is distributed by the National Optical Astronomical 
Observatory, which are operated by the Association of Universities 
for Research in Astronomy, under contract with the National Science 
Foundation. 



OB 


UT date 


UT„„ f 


exptime 


airmass 


seeing 




(Y-M-D) 


(h:m:s) 


(s) 




(arcsec) 


A 


2010-12-10 


06:26:33.886 


2345 


1.085 


0.51 


B 


2010-12-10 


07:19:18.192 


2450 


1.064 


0.48 


C 


2010-12-13 


05:07:19.888 


2450 


1.205 


0.66 


D 


2011-01-26 


04:14:19.507 


2450 


1.064 


0.55 


E 


2011-02-27 


02:59:58.770 


2450 


1.096 


1.26 


F 


2011-03-06 


02:35:28.955 


2345 


1.100 


1.30 
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Fig. 4. Histogram of all RVs for the FLAMES spectra (red, 
open histogram) and for the Besangon model (Rob in et al 



20031) in the direction of Be 39, on the same area of the 



Kassis et al. I (1 19971) data and normalised to the number of 



stars (filled, grey histogram). 
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Fig. 5. Six spectra obtained for star #2976, a binary near the 
base of the RGB. The right axis indicates the time since the 
first observation. 



3. Analysis 

3.1. Atmospheric parameters 

For this cluster we cannot apply the same proc edure as we 



adopted for the GC survey dCarretta et al. . l2Q09allbh . First, we 



saw that the photometry alone cannot guarantee atmospheric 
parameters with the required precision (at least, not for all 
stars). Second, we cannot use the IR K magnitude as we have 
previously done for other clusters to determine the tempera- 
ture from the V - K colour. Our sample does not only include 
bright stars for which the 2MASS photometry would be precise 



enough, but also faint stars for which it would be unreliable. We 
could potentially infer false inhomogeneities were we to use 
this colour. On the other hand, B, V, and / photometry accu- 
rate enough (apart from the mismatches discussed in Sec. I2.U 
is available for all targets. We therefore opted for temperatures 
derived from B - V and V - I as input values. 

We used the reddening and distance modulus derived 
above, adopting E(V-I) = l.3E(B-V),A v = 3AE(B-V), and 
4.75 as bolometric magn itude of the Sun. We u s ed th e colour- 
temperature relations of iRamfrez & Melend ez (2005) to find 
r e ff, and averaged the values obtained from B - V and V - I. 
The gravi ty was derived as in previous papers following the re- 
lations o f Alonso et al.l (ll999) for the bolometric correction and 
a mass of 1.15 M , as deduced from the isochrones. Average 
photometric T e g and logg values were computed, weighting 
more those obtained from the original colours than the "cor- 
rected" ones (see Sect. 12. 111. 



These r e fif's and logo's were used to derive first-pass iron 
abundances for all stars in our sample. The microturbulent ve- 
lociti es v t were obtained from the relation v, = -0.322 logg + 
2.22 ( Gratton et al. 19961) . Adopting this relation, only very 
weak trends in the relation between abun dances from F e neu- 
tral lines and expected line strength (see Magain, 1984) were 
apparent, with negligible effects on the abundances. We then 
computed the T e s and log g values required to satisfy the exci- 
tation and ionisation equilibria (we recall, however, that there 
are only three useful Fe n lines in the GIRAFFE HR13 wave- 
length range). The final adopted parameters are the average 
of the spectroscopic and photometric ones. Metallicities for 
all stars w ere then obtain ed with model atmospheres interpo- 
lated in the Kurucz ill 9931) grid of model atmospheres (without 
convective overshooting) with atmospheric parameters whose 
abundance matches that derived from Fe i lines. 

Temperature, gravity, microturbulent velocity, and Fe i and 
n are given in Tab. ?? for the UVES and Giraffe bright sam- 
ple (RGB and RC stars). The run of iron abundances with T e g 
and log g are shown in the upper panels of Fig. [6] These abun- 
dances do not display any significant trends with T e $ (left-hand 
panels) or log g (right-hand panels). The weak trend that is pos- 
sibly presented by the lower gravity stars in the [Fe/H] i-log g 
plane is negligible in the context of the present paper and could 
be caused by a variety of sources, such as blends, continuum 
tracing, or inadequacy of the model atmospheres - all more no- 
ticeable at cooler temperatures. We note, however, that there 
seems to be a small offset between the mean iron abundances 
of UVES and GIRAFFE stars. This is most probably due to the 
combination of higher resolution (hence accuracy in equivalent 
width, EW, measurement, see next section) and higher S/N for 
the UVES spectra. It is not worrisome, but in the next sections 
we will keep separate abundances from the two samples. 

The values for the SGB and TO stars have larger errors 
and larger scatter because of the low S/N combined with the 
warm temperatures and the larger errors in the photometry. 
Furthermore, we could not derive O abundances for the faint 
sample: the stars are too warm to measure the [O i] line, but 
not warm enough for the triplet at 6155-6158 A. While we pro- 
ceeded with a preliminary analysis also for the faint sample, 
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Fig. 6. Run of Fe i, Fe n, Na, and O with r e ff and \ogg for stars on the RGB and RC. Filled blue symbols are for UVES, open 
orange squares for GIRAFFE targets. 



we eventually decided to discard these stars and keep only the 
bright sample, for which our conclusions are sounder. 



3.2. Equivalent widths, sensitivity, and errors 



We m easured EWs wi th the software pack age Rosa ( GrattonL 



1988), as described in iBragaglia et alJ (120011) . We adopted the 
same automatic procedure to define the local continuum around 
each line as in previous papers, but we did not correct the 
GIRAFFE EWs to the UVES system because we did not have 
any star in common. Small differences between the two sets of 
EWs related to the different resolution of the spectra may in- 
deed exist, and explain the small differences in the abundances. 
The EWs of the member stars are made available only through 
the CDS. 



To evaluate the sensitivity of the derived abundances to 
the adopted atmospheric parameters we repeated our abun- 
dance analysis by changing only one parameter each time. The 
amount of the variations in the atmospheric parameters and the 
resulting response in abundance changes of Fe, O, Na, and all 
elements measured (i.e., the sensitivities) are shown in Tab. [4] 
In the upper part of the same table we also give the error con- 
tributed by each parameter. The column labelled "total error" 
gives the total error expected from uncertainties in the atmo- 
spheric parameters and in the EWs, to be compared to the 
observed scatter in order to understand if we have a statisti- 
cally significant variation in any of the elements. Note that for 
some elements (including Ca, V, and Ba) the observed scatter is 
larger than the observational errors. However, for various rea- 
sons errors may have been underestimated for these elements: 
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Table 4. Sensitivity to errors 



Element 




logg 


[A/H] 


v, 


EW 


Total 


Observed 




(K) 


(dex) 


(dex) 


km s _1 




error 


rms 






A parameter 












100 


0.30 


0.10 


0.2 












Error 












30 


0.03 


0.04 


0.1 








Fe i 


0.071 


0.014 


0.011 


-0.093 


0.013 


0.053 


0.049 


Fe ii 


-0.164 


0.046 


0.026 


0.043 


0.036 


0.066 


0.067 


Oi 


-0.058 


0.114 


0.027 


0.087 


0.085 


0.098 


0.119 


Nai 


0.004 


-0.025 


-0.015 


0.059 


0.052 


0.060 


0.050 


Mgi 


-0.032 


-0.008 


-0.007 


0.071 


0.053 


0.065 


0.066 


All 


-0.004 


-0.019 


-0.013 


0.069 


0.054 


0.064 


0.025 


Si i 


-0.101 


0.033 


0.010 


0.061 


0.047 


0.064 


0.065 


Cai 


0.035 


-0.063 


-0.014 


-0.008 


0.019 


0.024 


0.093 


Sc ii 


-0.087 


0.118 


0.025 


0.028 


0.053 


0.063 


0.060 


Til 


0.079 


-0.020 


-0.021 


0.003 


0.036 


0.044 


0.053 


Tin 


-0.103 


0.120 


0.023 


0.074 


0.106 


0.117 


0.052 


Vi 


0.096 


-0.012 


-0.019 


0.030 


0.023 


0.041 


0.129 


Cri 


0.091 


-0.008 


-0.024 


0.011 


0.087 


0.092 


0.150 


Mm 


0.031 


-0.052 


-0.001 


-0.039 


0.044 


0.050 


0.021 


Co i 


-0.027 


0.036 


0.008 


0.062 


0.090 


0.096 




Nil 


-0.037 


0.040 


0.007 


0.033 


0.019 


0.028 


0.045 


Ban 


-0.046 


0.063 


0.033 


-0.084 


0.046 


0.066 


0.120 



lines of Ca and Ba are strong and saturated, and those of V are 
affected by strong HFS effects. 

Errors in our measures are the following. 

J e ff '■ The errors in T e $ were derived by examining the rms of 
residuals along a quadratic log g - T e ff relation for RGB stars 
(after eliminating RC stars). This residual is 29 K, so we as- 
sume an error of ±30 K. 

logg: The error in logg is dominated by the errors in T e ff (30 
K, yielding an error of 0.006 dex) and in photometry (given 
all the uncertainties noted in Sec l2.11 we assumed 0.05 mag, 
yielding an error of 0.02 dex). This yields a total error of ±0.03 
dex. 

[A/H]: The error is not larger than the observed star-to-star 
scatter, i.e., ±0.04 dex. 

v,: This is the dominant source of errors for iron. It was es- 
timated by considering the star-to-star scatter of the trends of 
abundances with expected line strength: the error is 0.09 km/s. 
We assumed an error of ±0.10 km/s. 

4. Results 

Abundances of all elements were derived from EWs (and from 
synthetic spectra for O). Hyperfine splitting corrections were 
applied for Sc, V, and Mn. We used the same lines and line 
parameters as in the p apers of the FLAM ES GC surv e y (e.g ., 



The average values for all elements measured are presented 
in Tab. |5J separately for the UVES and GIRAFFE samples, 
given the diversity of resolution, wavelength coverage, and 
number of lines available in the two cases. For UVES the aver- 
age was derived from seven stars; for GIRAFFE it is typically 
derived from 21 stars. As stated above, only the brighter sam- 
ple was considered. Furthermore, star 1 1 was always excluded 
because is more that 3<x away from the cluster average veloc- 
ity and its approximately solar metallicity is clearly higher than 
that of the other stars of Be 39. 

Be 39 clearly is a homogeneous cluster: the star-to-star 
scatter is always explained by the errors expected from un- 
certainties in atmospheric parameters and EWs (see Tabs. [4] 
and|5]l. In particular, we find a mean metallicity [Fe/H]=-0.23 
(rms=0.04) from UVES and -0. 18 (rms=0.06) from G IRAFFE 
spectra. This agrees very well with Friel et al. ( 2010h (see also 
Sect. 14.2b . It is also in accord with metallicity expectations 
for clusters at Be 39's Galactocentric distance; the cluster lies 
at an Roc of about 11 kpc, where OCs generally have sub- 
solar metallicity, as shown in Fig. [7] Interestingly, its Rqc is 
around the transition between inner and outer disc, where the 
metallicity gradient begins to flatten. This is a region where 
more clusters should be obser ved to better constrain the mod - 
els of chemical evolution (e.g. jMagrini et al. , 2009; Friel et al 



irveyj 

Carretta et all l2009bl) . originally from iGratton et ail d2003b . 



Reference solar abundances are explicitly given in Tab. [5] 



l2010h . For the figure we selected only clusters older than 100 
Myr and with metallicity derived by high-resolution spectra; 
the sample contains 79 old OCs. The picture is quite clear; 
we caution, however, that the parameters (age, distance, metal- 
licity) have not been put on a homogeneous scale since this 
would be a great effor t, completely o utside the goal of the 
present paper (see e.g., lYong et aL 20121 for a recent work on 
the abundances and trends in the outer disc.) We also confirm 
Friel et al.'s derivation of a mild overabundance of a-elements 
([ar/Fe] =0.1.) 

Elements Al, Ti n, and Mn could only be measured in the 
UVES spectra. It is unfortunate that we could not measure Al 
also for the GIRAFFE sample, since it is involved in the light 
elements (anti-)correlations a nd may show also v ery strong 
star-to-star variations (see, e.g JCarretta et al. . 12012 . for the GC 
NGC 6752). But the UVES results suggest that this is not the 
case for Be 39 (see Tab.|5]l. 

Barium is the only neutron-capture element we measured. 
Its abundance was derived from three lines in the UVES spectra 
at 5853.69, 6141.75, and 6496.91 A. These lines yield abun- 
dances that agree well, with a typical rms of 0.12 dex for 
any given star. The star-to-star variations (0.08 dex) are due 
to the uncertainties in atmospheric parameters, especially in 
T e g and v f . We have only a single Ba line (at 6141.75 A) in 
the GIRAFFE spectra. There is an offset between the Ba abun- 
dances measured in the two types of spectra, which can be ex- 
plained by differences in the EWs, however. The EWs of the 
6141 A are larger for GIRAFFE than for UVES stars of simi- 
lar atmospheric parameters, most probably because of an unre- 
solved blend with Ni i and Si i in the lower resolution spectra. 
The contaminant-line contributions to the GIRAFFE Ba lines 
seems confirmed by estimates of their fractions to the total EW 
of the blended Ba n feature. 
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Fig. 7. Distribution of [Fe/H] with Galactocentric distance for 
OCs, colour-coded according to age (light blue triangles : age 
0.1-1 Gyr; orange squares: age 1-4 Gyr; red circles: age older 
than 4 Gyr). [Fe/H] and errorbars are taken from the original 
literature sources. Be 39 is indicated with a filled red symbol. 



Given the strength of the Ba n lines, which are always sat- 
urated, Ba is a difficult element to measure. However, three 
stars of the faint sample (1082, 1357, and 2170) show a very 
high Ba abundance; all three have RVs compatible with cluster 
membership. The first is the most convincing case suggestive 
of an s-process-enriched "Barium star". Not only is its Ba line 
much stronger than in stars of similar atmospheric parameters 
and normal Ba abundance (see Fig. [8]), it also shows a clear La 
ii line at 6320.4 A, not detected in other, similar objects. Star 
1082 is brighter than expected from a single star at the TO and 
could be a blue straggler; the high Ba could be due to mass 
transfer in a binary system. Since stars on the main sequence 
have only a tiny convective envelope, it does not need to have 
accreted a lot of mass to explain the higher Ba abundance; this 
would be possible even in a wide system and this would explain 
why its RV is not discrepant from the cluster' average value. 

4.1. Na and O abundances 

The main goal of our study was determining the Na and O 
abundances in Be 39 members to see whether they show an in- 
trinsic star-to-star scatter or even an anti-correlation, the main 
spectroscopic signature of multiple populations in clusters. 

We measured oxygen abundances from the [O i] line us- 
ing spectrum synthesis. W e considered the Ni i contamina- 



tion of the line (see e.g., ICarretta et aU I2007I) . as well as 



the coupling of O and C expected in cool stars. The disso- 
ciation equilibrium was computed assuming [C/Fe]=-0.2 and 
[N/ Fe]=0.5, which are typi cal values for red giants in old OCs 



(see iMikolaiti s et aUl2012l) . To estimate the sensitivity of our 



O abundances to this assumption, we note that they would be 
higher by ~ 0.05 dex for RC stars and by 0. 1 1 dex for stars at 
the tip of the RGB had we instead assumed [C/Fe]=0. 

Only the Na doublet at 6154-6160 A was used to mea- 
sure sodium abundances since the Na D and 5686-5690 lines 



Fig. 8. Comparison of the Ba n 6141.7 A line (indicated by an 
arrow) in the Ba-star 1082 and in the star 1468, of very close 
atmospheric parameters, as confirmed by the similarity of the 
other lines seen here (Ca, Ni, two Fe i, Fe i+Fe n, Fe n, Fe i, 
Na, and Si, from left to right and indicated by grey lines). The 
spectra are shifted for clarity. 

appear only in the UVES spectra and are too strong at this 
metallicity. These abund ances were correc ted for NLTE ef- 
fects using the relation bv lLind et al.1d201#l We interpolated 
their corrections with the quadratic relation A^ue - Aue = 
-0.088 + 0.00076394 x EW - 0.0000144444 x EW 2 , valid for 
EW< 150 mA (that is, all Na lines in our stars). The [Na/Fe] ra- 
tios (see Tab. ??) were then computed assuming [Fe/H] =-0.23 
for UVES stars and -0.18 for GIRAFFE stars. 

Neither O nor Na abundances show any significant trend 
with T e ff or log g (see lower panels of Fig . |6j> . As apparent from 
Fig- and from the very small scatters associated both to the 
O and Na values (Tab. |5}, Be 39 is very uniform in these ele- 
ments as well. The star-to-star scatter is completely explained 
by the associated internal errors; the only effect mimicking a re- 
lation between O and Na is the small offset between the UVES 
and GIRAFFE samples. Even with all the caveats already pre- 
sented, this is a solid result: Be 39 is a normal OC, without 
any signature of the Na- O anti-correlation, w hich is the typical 
"tag" for a genuine GC dCarretta et al.Ll201db . 



4.2. Comparison with lFriel ef~a/1 \2010\) 

Friel etalJ d2010t) observed two RC and two bright RGB stars 
using the echelle spectrograph on the KPNO 4m telescope. We 
give the data for the three objects in common with our study in 
Tab. |6j two of them were observed with UVES fibres, one with 
GIRAFFE. The resolution of their spectra is R ~ 28000, inter- 
mediate between our GIRAFFE and UVES spectra, the wave- 
length coverage is AA = 4800 - 8100 A, and the S/N is always 



3 We generally use the iGratton et al] {l999) corrections in our 
works on GCs; however, g iven the much hi gher metallicity of Be 39, 
we preferred to adopt the iLind et al] d201 lh method. This does not 
mean that we need to revise all our Na abundances, since the older 
NLTE corrections were computed for the metallicity regime of GCs. 
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0.4 - 



0.2 



-0.2 



-0.4 




Table 6. Comparison with iFriel et al.1 d20ld) for two UVES 
stars and one GIRAFFE star. 



Fig. 9. Na and O distribution for the UVES (filled blue points) 
and the bright GIRAFFE (open orange squares). The rms in 
O and Na are shown by errorbars of different colour. The ex- 
pected errors on O and Na abundances (see Tab. |4|i are shown 
as black errorbars in the left lower corner of the plot; they are 
comparable to the observed rms. 

Table 5. Average abundances for Be 39. 



element 


Sun 


abu 


rms 


abu 


rms 


Notes 






UVES 


Giraffe 




[Fe/H] i 


7.54 


-0.23 


0.04 


-0.18 


0.06 




[Fe/H] ii 


7.49 


-0.19 


0.05 


-0.21 


0.07 




[O/Fe] i 


8.76 


0.03 


0.07 


-0.08 


0.11 




[Na/Fe] i 


6.30 


-0.02 


0.04 


-0.01 


0.03 


NLTE 


[Mg/Fe] i 


7.43 


0.17 


0.02 


0.24 


0.05 




[Al/Fe] i 


6.40 


0.02 


0.03 








[Si/Fe] i 


7.53 


0.12 


0.03 


0.14 


0.05 




[Ca/Fe] i 


6.27 


0.07 


0.05 


0.17 


0.08 




[Sc/Fe] n 


3.13 


0.00 


0.05 


-0.03 


0.06 


HFS 


[Ti/Fe] i 


5.00 


0.03 


0.03 


0.01 


0.06 




[Ti/Fe] ii 


5.07 


0.12 


0.05 








[V/Fe] i 


3.97 


-0.08 


0.04 


-0.16 


0.08 


HFS 


[Cr/Fe] i 


5.67 


-0.15 


0.05 


-0.13 


0.08 




[Mn/Fe] i 


5.34 


0.00 


0.02 






HFS 


[Ni/Fe] i 


6.28 


0.04 


0.04 


0.07 


0.03 




[Ba/Fe] n 


2.18 


0.14 


0.08 


0.30 


0.09 





high, between 70 and 115. Friel an d coworkers noted that one 
star (WEBDA ID 2619, or 2130 for lKassis etafl 19971 adopted 
as our ID) is pro bably a binary, comb ining their data and the 
RV measured by Frinchabov et al .(2006), differing by about 10 
km s . We do not see any strong indication of an RV change 
in our data, but the average RV is about 1.5 km s _1 lower than 
theirs. A similar behaviour is shown by star 1435/1139, again 
with our RV lower by about 6 km s _1 but no strong signature of 



Star 


T eS log g v, 


RV 


ID 


T eff log g v, 


RV 


U/G 




Friel et al. 2010 






This paper 






2055 


4450 1.8 1.5 


57.5 


1657 


4399 2.10 1.54 


57.5 


U 


2619 


4750 2.2 1.5 


60.0 


2130 


4740 2.43 1.44 


58.5 


u 


1435 


4750 2.2 1.5 


62.8 


1139 


4659 2.46 1.43 


56.9 


G 



Element 


abu 


rms 


abu 


rms 


Friel+10 




UVES 


Giraffe 




[Fe/H] i 


-0.23 


0.04 


-0.18 


0.06 


-0.21 


[O/Fe] i 


0.03 


0.07 


-0.12 


0.11 


0.02 


[Na/Fe] i 


-0.02 


0.04 


-0.01 


0.03 


0.09 


[Mg/Fe] i 


0.20 


0.02 


0.22 


0.05 


0.19 


[Al/Fe] i 


0.05 


0.03 






0.20 


[Si/Fe] i 


0.14 


0.05 


0.14 


0.05 


0.20 


[Ca/Fe] i 


0.10 


0.05 


0.15 


0.08 


0.01 


[Ti/Fe] i 


0.06 


0.03 


-0.01 


0.06 


0.03 


[Ti/Fe] ii 


0.15 


0.05 






-0.09 


[Cr/Fe] i 


-0.13 


0.05 


-0.13 


0.08 


0.06 


[Ni/Fe] i 


0.07 


0.04 


0.05 


0.03 


-0.06 



RV variation (only one discrepant value in our six measures). 
We therefore define those two stars as probable binar ies. 



Tab. [6] displays a direct comparison of the IFriel et al 
d2010h atmospheric parameters and ours; they gen- 
erally agree well. We also compared our EWs with 
theirs; the agreement is excellent, with our EWs slightly 
smaller. We find EWuves = 0.992(+0.019) x EW Frie i - 3.86 
mA (with rms=0.64 mA, over 102 lines) and 
EWgiraffe = 1 .030(±0.046) x EW Frie i -1.1 mA (with 
rms=0.76 mA, over 26 lines). 

A detailed comparison of the two analyses is beyond the 
scope of our paper. We only note that the derived abundances 
agree well (see Tab. [6]), especially for iron and the a-elements. 
We also recall that we corrected the Na abundances for NLTE, 
while Friel and collaborators did not, since these corrections 
are small for such high-metallicity, warm stars. Finally, we do 
not provide here any comparison with the properties of other 
clusters, or discuss the radial abundance distribution. This is 
outside the main goal of our stu dy; moreover, th e abundances 
are so similar to those derived by Friel et al. ( 2010b that all con- 
siderations and conclusions of these authors apply. 



5. Summary and conclusions 

We presented the abundance analysis of the high mass, old OC 
Be 39. The high-resolution spectra of more than 100 stars, ob- 
tained with FLAMES @VLT, permitted us to isolate about 70 
cluster members. With the advent of surveys targeting large 
numbers of stars in large and significant sets of clusters (e.g., 
the Gaia-ESO 'Surve y, recently started on FLAMES @VLT, see 
Gilmore et alJl2012h . comparable samples will become more 
common, but this is presently the largest sample of stars in this 
particular cluster and one of the largest in general. The mem- 
bership information was used to determine updated cluster pa- 
rameters (age, reddening, and distance). About one half of the 
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stars were analysed to determine metallicity and detailed ele- 
mental abundances. We confirmed the slightly sub-solar metal- 
licity of Be 39 ([Fe/H]=-0.2) and the very small scatter around 
the average abundances for all elements measured (Fe, O, Na, 
Al, a-, iron-peak, and heavy elements). This is interesting also 
in the context of cluster and disc formation mechanisms and 
disc chemical evolution. 

The observations were optimised to look in particular for 
possible star-to-star variations in O and Na, which are indica- 
tive of the existence of multiple populations in this cluster, in 
analogy to what is is found for the higher mass, older GCs. No 
such scatter or anti-correlation was found; and we conclude that 
Be 39 is a normal, homogeneous, single-population cluster. 
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